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ABSTRACT 

A packed bed reactor model for the sulphation of active sodium carbonate has been 
developed to simulate sulphur dioxide breakthrough curves. Breakthrough cur\'es 
obtained from the formulated model were then validated by comparing them with those 
available in literature. Kinetic data and other physical parameters available in the 
literature were used to predict the breakthrough curves for the above mentioned system 
using the formulated model. Simulated breakthrough curves corroborated fairly well with 
the experimental results reported in the literature. 
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NOMENCLATURE 


A 

Ag 

Ap 

B.D 

b.c.d 

Ca 

Cao 

De 

De(0) 

De(X) 

Fp,fp 

k 

ro 

t 

t*,T’ 

T 

U 

Vg 

Vp 

w 

y 

* 

y 


Reactant gas (SO2) 

Surface area of grains of solid, cm’.'sec 
Area of a particle, cm" 

Solid reactant and product respectively 
Stoichiometric coefficients 

Molar concentration of SO2 at the outlet, g.moles/cm 

Molar concentration of SO2 at the inlet, g.moles/cm' 

Effective diffusivitv'. cm", sec 

Initial effective diffusivit\% cm" 'sec 

Effective diffusi\dt\' at conversion x. cm'/sec 

Shape factor for grains of solid 

Shape factor for particle 

Intrinsic reaction rate, cm/sec 

Grain radius, cm 

Initial radius of grain, cm 

time of reaction 

dimensionless reaction time 

tempreture, K 

superficial velocity, cm/'sec 

Volume of grains of solid reactant, cm^/gm 

Pore volume of the particle, cm^ 

Fractional volume of reactive solids 
Fractional conversion of solid 
Distance along the flow, cm 
Dimensionless distance 



Gr^it Symbols 


a A distributed parameter 

jS A function of x and a 

Y Grain size change factor 

Porosity of a particle 
Porosiw at x 

.4 

Sq Initial porosity' 

e,. Bed porosity' 

^ Ratio of ty and Tq 

a Characterstic reaction modulus 

Pg Molar density' of solid reactant 

¥ 


Dimensionless concentration 



1. INTRODUCTION 


The increased presence of sulphur dioxide in the atmosphere causes several adverse 
effects on human health, vegetation, material and visibility. Excessive emissions of 
sulphurdioxide lead to acidification of rain water, lakes and soils. The various methods 
for control of sulphur dioxide emission include burning less sulphur fuels, removing 
sulphur from fuel, substitution of another energy source, cleaning up of the flue gas for 
removal of sulphur dioxide as far as possible right at source and eventual dispersion 
through stacks. 

-Among various control methods, cleaning up of flue gases received maxmum attention in 
recent years. The increased interest in dry flue gas desulphurisation processes has led to 
the need for evaluating various sorbents which remove sulphur dioxide by chemical 
reaction, adsorption, catalytic reaction or a combination of processes. Sulphur dioxide 
removal by porous solid reactants was found to be limited by one or more of the 
following steps: Bulk diffusion of sulphur dioxide to the particle, diffusion of sulphur 
dioxide through pores, diffusion of sulphur dioxide through product layer, filling up of 
small pores causing decrease in reactive area, built up of product layer at entrance of 
pores causing pore closure, less surface area due to sintering and reduced kinetic reaction 
rates at low tempretures. 

Evaluation of sorbents under various process conditions may be facilitated if their 
performence can be predicted theoretically, using an appropriate model. In the present 
study an attempt has been made to formulate a mathematical model to predict 
breakthrough curves for flue gas desulphurisation under dry condition in a packed bed 
reactor using active sodium carbonate as a sorbent. Before modelling an existing physical 
system, it is essential to examine it and then decide parameters which can be used to form 



mathematical expression representing the physical problem. Breakthtough curv'es 
obtained from the formulated model can then be validated by comparing them with the 
actual results obtained from the physical system which was modelled. Once the simulated 
breakthrough curves are validated they can be used to predict the performence of suitable 
sorbents in a packed bed reactor under simmilar conditions. The same can then be used to 
design packed bed reactors for the removal of sulphur dioxide using suitable sorbents. 



2. BACKGROUND INFORMATION 


2.1 Sulphur Dioxide Control Technologies 

There are four possible approaches that may be used individually or in combination to 
control sulphur dioxide emission from fossil fuel combustion. They include fuel 
switching, fuel cleaning and fuel conversion, desulphurisation of flue gases, dispersion 
from tall stacks. Several combustion power facilities around the world ha\'e switched 
over to low sulphur fuels which is found to be fastest and easiest means of complying 
with sulphur dioxide emission regulations. Sulphur removal from fuels before 
combustion represents a more likely alternative if technology for removing sulphur from 
fuels is available for economic implementation. Sulphur is present in coal in varying 
amount in pyritic and organic forms. Tall stacks are obvious means of reducing ambient 
ground level concentration of sulphur dioxide at cost significantly less than those for flue 
gas scrubbing but total amount of emission remains the same which may adversely affect 
the environment in many ways. 

An optimal combination of flue gas desulphurisation for the control of sulphur dioxide at 
source and subsequent dispersion with the help of suitable stack might be described as 
the approach for achieving pollution reduction relatively independent of uncertainties of 
meteorology and fuel quality. Flue gas desulphurisation is presently the most commonly 
used technology to comply with sulphur dioxide emission requirements. In fact dr>' flue 
gas desulphurisation will be a promising technology compared to conventional wet 
scrubbing in view of the requirement of dispersion as stated above. 

The flue gas desulphurisation may be grouped according to two classifications 

1 . Throw away or regenerative 

2. Dry or wet 
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Throw away process is that in which waste product is discarded, as a result fresh 
chemical must be continually added. In regenerative process as the name implies, the 
chemistry is such that the removal agent can be continually regenerated in closed loop 
system . Wet or dry processes are differentiated simply by whether the product formed is 
dry or wet. The removal system typically involves the use of adsorption, absorption or 
catai>l:ic processes. 

Dr>' desulphurisation of flue gases may be achieved by passing flue gas thorugh the 
packed bed containing reactive sorbent dispersed in an inert media. Desulphurisation 
takes place due to reaction of sorbent with sulphur dioxide, knowledge of gas-solid 
reaction is thus essential in order to understand the phenomena of sulphur dioxide 
removal using reactive sorbents. 

2.2 Elements of Gas Solid Reaction Systems Involving Single Particle 

The smallest representative unit of gas-solid reaction system is the interaction of single 
particle with a moving gas stream. Let us consider a gas-solid reaction of the type 

A(g) + bB(s)=cC(g)+dD(s) 11 

Where b,c and d are the stoichiometric coefficients. Over all reaction process may invok e 
the following individual steps: 

1 . Gas phase mass transfer of the gaseous reactant from the bulk of the gas stream to the 
external surface of solid particle. 

2. Diffusion of gaseous reactant through the pores of the solid matrix, w'hich could 
consists of a mixture of solid reactants and products. 

3. Adsorption of the gaseous reactant on the surface of solid matrix. 

4. Desorption of the gaseous product from the surface of solid matrix. 

5. Diffusion of gaseous reactant product through the pores of solid matrix. 



6. Gas phase mass transfer of the gaseous product from the external surface of the solid 
to the bulk of the gas stream. 

Heat transfer: 

For exothermic or endothermic reactions the diffusional and reaction steps will be 
accompanied by 

1 . Convective (and possibily radiative ) heat transfer between the gas stream and surface 
of solid particle, and 

2. Conduction heat transfer within the solid reactant-product matrix. 

Structural changes: 

The reaction and heat transfer processes could lead to structural changes, such as 
sintering or changes in pore structure, which in turn could have a marked effect on the 
over all reaction rate. 

2.3 Grain Model 

In many gas-solid reactions encountered in chemical and metallurgical processes, solid 
products are formed. Such reactions may be described by the general equation 

A(g) bB(s) = cC(g) + dD(s) 2.2 

The overall volume of solid may increase or decrease depending on the relative densit\- of 
solid product compared with that of solid reactant. In most cases however the change is 
rather small so that the overall size may be regarded constant. 

When solid reactant is non porous there is a sharp boundary between the reacted and 
unreacted core. In general case of porous solid there is a gradual change in the degree of 
conversion throughout the particle. The external layer will be completely reacted after a 
certain time, and the thickness of completely reacted layer will increase towards the 
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interior of the particle. Under these conditions in contrast to non porous solids the 
reaction within the partially reacted zone occurs simultaneously with diffusion of fluid in 
this zone. 

When chemical reaction presents the major resistance to the overall progress of reaction, 
the concentration of fluid reactant will be constant everywhere and the reaction will occur 
uniformly throughout the volume of the solid. If on the other hand, pore diffusion 
presents the major resistance, the reaction will occur in narrow boundary between the 
unreacted and completely reacted zones. In the intermidiate region where the resistances 
presented by chemical reaction and pore diffusion are of comparable magnitude, both 
processes must be considered simultaneously. 

2.3.1 Assumptions Made in the Formulation of Grain Model 

A particle is usually made up of fine grains with or without some binding agents. The 
overall shape may be approximated by that of a slab, a long cylinder or a flat plate, made 
up of individual grains of equal size, which could again be spheres, long cylinders or flat 
plates. In cosidering such a particle to be made up of large number of discrete grains, we 
are not far from the truth, although in genial the assumptions of uniform size of grain and 
its regular shape are not accurate. 

The assumptions made in the formulation of grain model are discussed by Shon and 
Szekely ( 1972) and will be merely listed here. 

1 . The pseudosteady state approximation is appropriate for describing the concentration 
of reactant within the particle. 

2. The resistance due to external mass transfer is negligible. 

3. Diffusion ■within the particle is either equimolar counterdiffusion or is at low 
concentration of diffusing species, and the effective diffusivities of gaseous reactant 
and product are equal and uniform throughout the particle. 
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4. The system is isothennal. 

5 . The solid structure is macroscopically uniform and unchanged by reaction. 

6. Diffusion of the gaseous reactants through the product layer of individual grain is not 
rate limiting. 

7. Viscous flow contribution to the mass transfer within the particle is negligible. 

2.3.2 Chemical Reaction Controlled Model 

Chemical kinetics control the process when the rate of diffusion through the interstices 
among the grain presents the negligible resistance to the progress of reaction. The 
concentration of gaseous reactant is uniform throughout the solid. The grains are assumed 
as non porous and within each grain the reaction front retains its original geometric shape 
as reaction proceeds. When the product layer has any degree of porosity, for the n'" order 
reaction it can be written as ( Shon and Szekely, 1972 ). 

T* = (l -( 1-X)‘^^^) 2.3 

* 

Where T , X, fr are dimensionless time, fractional conversion and grain shape factor 
respectively. The grain shape factor has value of 1 , 2 and 3 for flat plates, long cylinders 
and spheres respectively. 

2.3.3 Reaction Controlled by Diffusion through the Product Layer 

In case of pore diffusion being rate limiting step, reaction occurs in a narrow zone 
separating the unreacted core and completely reacted layer, where the reactant 
concentration becomes very small. Under this situation it can be written as ( Shon and 
Szekely, 1972 ) 


T* = o-^ A(X) 


Where, 

A(X) = X^ Forfp = l 


2.4 


2.5 
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= X + ( 1 - X ) In ( 1 - X ) For ^ = 2 2.6 

= ( 1 - 3 ( 1 - X ) + 2 ( 1 - X ) ) For fp = 3 2.7 

and cr is the reaction modulus 

2.3.4 Reaction Controlled by Simultaneously Chemical Reaction and Pore 
Diffusion 

When chemical reaction and pore diffusion control the overall rate of reaction then 
approximate relationship has been proposed by Shon and Szekely ( 1972 ). 


T* = ( 1 - ( 1 - X ) + a - A(X) 


Where ct is a charecterstic parameter called reaction modulus. Its value is given as 


a = 


Ze 

A„ 


, 2D, ^ 


2.9 


An estimation of reaction modulus provides information as to one needs data on chemical 
kinetics or pore diffusion or both, in order to predict the relationship between conversion 
and time. If <j is much smaller than unity, one needs kinetic parameter only, if it is much 
larger than unity , than only diffusion parameter is required. If g is around unity, both 
kinetic and diffusion parameters are required. 

2.3.5 Diffusion Controlled Model with Changing Porosity 

In some system, the reaction of porous solid with a gaseous reactant often seriously 
affects the macroporosity of particles due to shrinking or swelling due to density 
differences of solid reactant and product. Accounting this swelling or shrinking by a 
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characterstic factor y ( Ramachandran and Smith, 1977 ) the pore diffusion model can 
be modified as 

r=cT/A(X) 2.10 


Where a^has been presented for varying diffusivity. Dj(x), which can be given as 
(Ghosh. 1987); 


De(x) = D,(o) 





2.11 


Where De(o) is the original diffusivity and is given as 

e.= l-0-eo)[y-(l-Y)0-X)] 2.12 

and 

Where b and d are the stoichiometric coefficients and pg and pj are the molar densities 
of solid reactant and product. Ify is greater than 1 swelling takes place and if less thanl 
shrinkage develops. 


2.4 Kinetics of Sulphation of Active Sodium Carbonate 


The active sodium carbonate is obtained by thermally decomposing pure sodium 
carbonate, the reactant becomes highly porous. Erdos et al. (1989 ) reported that active 
sodium carbonate prepared by thermal decomposition of sodium carbonate is more 
reactive than anhydrous sodium carbonate prepared by dehydration of hydrated 
carbonates. 
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Thermal decomposition of sodium bicarbonate above 115°Cis reasonably fast at 
atmospheric pressure yielding chemically anhydrous sodium carbonate according to the 
equation 


2NaHC03 = Na2C03 + CO2 ^ H2O 


2.14 


Ff= 6.79 K Cal 

The tempreture at which decomposition occurs will depend on the partial pressure of 
carbon dioxide and water vapour present, which is well below equillibrium values at 
stack gas conditions. The sodium carbonate formed is known to have high surface area. 
The tempreture at which sodium bicarbonate decomposes has been shown to influence 
the reactivity of sodium carbonate with sulphur dioxide (Erdos et. ah). 

2.4.1 Reaction with Sulphur Dioxide 

The overall reaction of sodium bicarbonate with sulphur dioxide has following form 

2 NaHC 03 + SO2 + 1/2 O2 = Na2S04 + 2 CO 2 + 2H2O 2.15 

Ff= -69.66 K Cal 

This is generally accepted reaction scheme. It assumes that reaction is followed by the 
reaction between active sodium ' carbonate and sulphur dioxide. The active sodium 
carbonate reacts readily with sulphur dioxide yielding sodium sulphite as follows 

Na2C03 + SO2 =^Na2S03 + CO2 2. 16 

Ff= -13.00 K Cal 

and in the presence of oxygen it reacts with sulphur dioxide yielding sodium sulphate 
according to the reaction 



Na.COj + SO2 + I/2O2 = Na2S04 + CO2 


2.17 


Ff = -76.45 K Cal 

Since the active sodium carbonate was obtained by thermally decomposing pure sodium 
bicarbonate, the reactant becomes highly porous. For porous solid reactant, gas-solid 
reaction models based on grain structure of solid have been reported to be appropriate 
(Ghosh and Maiti, 1989 ). 

For conditions typical of the experiments for sulphation of active sodium carbonate, 
predicted reaction time, t, for complete conversion was calculated to be 85 seconds for 
pellet of 12.7 mm diameter using film diffusion model by Bhaskar. 1994. This predicted 
time was much less than the actual time required for same conversion as observed by 
Bhaskar, 1994, thus confirming little contribution by gas film diffusion for the overall 
rate. 

2.4.2 Effect of Tempreture on Sulphation 

In general sulphation of active sodium carbonate at high tempreture is preffered due to 
increased reaction rate at high tempretures. Increase in conversion due to increase in 

tempreture firom 250°Cto 300°Cis negligible (Bhaskar, 1994 ). This may be due to 
sintering effects of active sodium carbonate which offsets the increased reaction rate in 
that range 

2.4.3 Structural Examination of Active Sodium Carbonate 

The scanning electron micrograph of active sodium carbonate prepared by decomposition 
at tempreture 150°C, 250°Cand 350‘’Cshowed the formation of gas holes due to 
evolution of carbon dioxide and water vapour (Bhaskar, 1994). The particles of active 
sodium carbonate were smaller when decomposed at 150° C whereas they were relatively 
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larger when decomposed at 250 “C, whereas at 350 °C active sodium carbonate lost both 
its grain structure as well as gas holes as observ^ed by Bhaskar, 1994. This may be due to 
sintering that might have occured at high tempreture. 

When active sodium carbonate is allowed to react with sulphur dioxide it was obsen.ed 
by Bhaskar, 1994 that at tempreture of 250° C growth of sulphated layer occured as well 
as there was a decrease in size of particles and gas holes. The decreased size of pores 
may have contributed to pore diffusion controlling the progress of reaction. 

The existence of straight line relationship according to chemical reaction model does not 
necessarily mean that the progress is controlled by chemical reaction, because tliis 
straight line is very difficult to distinguish from the curve for pore diffusion in man\- 
cases. This leads to an incorrect assumption of chemical reaction control, ^vhereas the 
reaction is controlled by ash layer diffusion or was a mixed control. 

Sulphur mapping was uniformly dense throughout the area of pellet for radial position, it 
was less dense towards the core of the pellet also a definite interface in the concentration 
of sulphur was observed (Bhaskar, 1994 ).This means that the sulphation reaction is 
uniform and complete on the outer surface of the pellet and there was a narrow zone in 
which sulphur concentration drops to a small value. This supports the hypothesis that 
reaction is chemically controlled for some initial time during which the reaction is 
uniform throughout the particle leading to partial reaction of inner core. As the reaction 
proceeds, the the formation of sulphate decreases the pore diffusion which leads to the 
formation of seprate zones of completely reacted and partially reacted portion radialy 
inwards the particle. 



3. SCOPE OF THE WORK 


While going for modelling it is essential to examine the physical problem, which is to be 
modelled. Physical problem is then replaced by equivalent mathematical expression.Once 
the mathematical expression representing the physical problem is formulated, it can be 
solved by suitable technique. Results obtained from the mathematical model are then 
compared with the actual results of the physical problem. 

In the present study attempt has been made to 

1. Formulate a fixed bed reactor model for the sulphation of active sodium carbonate 
(freshly prepared by heating sodium bicarbonate in the reactor) dispersed in sand. 

2. Simulate breakthrough curves for the above mentioned system using the formulated 
model. 

3. Validate the simulated curves by comparing them with the experimental results 
available in the literature. 



4. PACKED BED REACTOR MODELLING 


Let us assume a fixed bed reactor filled with 
introduced from the bottom of the reactor. 


reactive porous solid. Pollutant gas is 


y+Ay 

y 

y = 0 


T 



T 


(t) 

O 

An elemental volume in the reactor is assumed whose thickness is 0’ ^nd cross sectional 
area is S. The pollutant gas is introduced in the reactor with a superficial velocity of U 
and is assumed to react with spherical porous reactive solid whose porosity’ is 
according to the grain model. The porous reactive solid is dispersed in a bed of porosity 
e,. .The fractional volume of the reactive solid of density in the reactor is w and the 
fractional conversion of the sorbent which vanes with reaction time is given as x. then the 
material balance can be written as follows 


Assuming U, S are constant, taking the limit <^t ->0 and rearranging the 

continuity equation takes the form 


(s, + e.)^c,s 


+ 


^SpgWX 


(l- e^)(l- €,.) 
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u 




5 y 


■ + 


/ c 


sc, , (l-e,) 


5 X 




4.2 


According to the grain model, each of a particle is reacting according to unreacted 
shrinking core model. The continuity equation given above may be expressed in terms of 


dimensionless parameters explained as follows containing the characterstic features of 
grain model. 


W 


C, 


C 




.■10 


P.F/. J 


and 


4.3 


T* 


- w 


V F 

g g 



)(1- ^v)T 


4.4 


writtmg the partial derivatives in terms of dimensionless variables, we have 


5 y 


5 ij/ k , 




(l- e^)(l- €,.) 


s~T 


S^( bkC^^Ag " 
5t\p,F^Aj 


C 


AO 


5 X _S X ^bkC ^f^Ag 
sl^FVlp.FgVg) 


4.5 


4.6 


4.7 


Substituting the derivatives into the continuity equation, final form may be obtained as 
follows 
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UC,^-w\ 




+ + e,) 




bkC^^ 

. J 


S y 
5 X 

Pb -.7 = o 


PbF/, 


c 

■Asr) 


4.8 


Eliminating the common term and rearranging the continuity equation we get 


5 y/ 5 w 5 X 

-_ + cr Z-+ =0 

by Sr 5 1 


4.9 


where 


a = 


b{^P + ^v)Cao 
w{l- €^(1- e,.)p, 


4.10 


Now X is a conversion factor which is given as 

X = (Initial volume of sosbent -volume of unreacted core ) / Initial volume of sorbent 


= 1 - 






where 




ii' / f 6 


4.11 


4.12 


^ = 


4.i; 


Therefore 


f7 = -^^)l7 


4.14 


4.15 


Now the local rate of reaction at the solid surface may be expressed as follows for a 


reaction, first order with respect to gaseous reactant at shrinking core surface (ty) 



5 

P B ~ bkC ^ 

o t 


17 


Therefore 


Sr^ _ -bkC, 

^ t Pb 

We can write 

^ ^ 

5 1' 5 t St’ 

Which can be written as 


Jr 
5 1' 


-bkC, 


Pb 


bk i A, 

pX^j 


c 


.^0 


4.16 


4.18 


Initial value of the radius of grain is given as 
therefore we can write 


FV 

g g 


S t’ 


= -hW 


4.20 


4.21 


5 X • 

or we can write — ^ = 3(4^)' = P (say) 


4.22 


An approximate value of fS with sufficient accuracy may be obtained from the following 
expression ( Shon and Szekely, 1972 ). 


P = 


1 


' +2cr'[(l-x)-l) 3-1 


Tl 


4.22 


The above expression is based on the assumption that pore diffusion and chemical 
reaction are simultaneously rate limiting. 
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In case of acti\ e sodium carbonate it was observ'ed that the theoretical results obtained 
from the model based on pore diffusion and chemical reaction being simultaneoush' rate 
limiting steps matched fairly well with the values of thermogravimetric analvsis for 
fractional conversion, x. by Bhaskar. 1994 (See^ppendix B ). 

The final equation is 


3 y/ 3 y/ . ^ 

^a--+/? = 0 
3 V dr 


4.24 


This mav be soh'ed with initial and boundaiv conditions 


y/{0j) = 1 = f[t'),y/[y\0^ = 0 = h(^y') 


4.25 


Assuming /? as constant the above differential equation may be solved by Laplace 


transformation (se^ppendix A). 


To evaluate values of fi^ctional conversion at^fier^nt dimensionless tim^Runge-Kutta 


/ ^ A 


method of fourth order was used (see^penc^ix E). 





5. RESULTS AND DISCUSSION 


5.1 Discussion on the Model Formulated 


The model formulated for the packed bed reactor assumes plug flow of gas. isothermal 
conditions in the bed of reactor, particles to be uniform and spherical in shape and 
negligible diffusive or dispersive mass transfer in the axial direction. Experimental results 
were taken from the literature (Bhaskar, 1994), in which fixed bed reactor was used 
containing Ig of active sodium carbonate dispersed in 40g of sand. An earlier smdy found 
that no sulphur dioxide removal occured with only sand in the reactor (Jozewicz and Co- 
workers. 1988). It wast therefore decided to use the term indicating the fractional volume 
of the reactive sorbent, w, present in the reactor while going for the packed bed reactor 
modelling. Although porosity may change with reaction time, in the present study it was 
assumed to be constant. The experimental data for fractional conversion, x. of active 
sodium carbonate with reaction time available in literature matched fairly well with 
combined pore diffusion and chemical reaction being rate limiting model rather than 
changing porosity model, see Appendix B. Therefore it was decided to keep pore 
diffusion and chemical reaction simultaneously rate limiting steps. Furthermore there is 
very little contribution by gas film diffusion for the overall rate in case of active sodium 
carbonate and sulphur dioxide system (Bhaskar, 1994), a constant reactant concentration 
in the bulk was assumed or more precisely, at the outer surface of the pellet. Which 
means dimensionless concentration outside the particle was assumed to be unity. In 
addition, for the values of reaction modulus, <J , more than 3 extent of reaction for the 
sorbent with respect to reaction time is nearly uniform throughout the length of the bed 
(Evans and Song, 1974).For active sodium carbonate and sulphur dioxide system at 
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tempretures 150°C,200°C and 250°C the value of a ranges between 3.45 to 3.67. 
Therefore it was decided to keep extent of reaction uniform throughout the bed. 

5.2 Simulation of Breakthrough Curves Using the Formulated Model 

Sulphur dioxide breakthrough curves for sulphation of active sodium carbonate dispersed 
in sand were simulated using the formulated packed bed reactor model. For the sulphation 
of active sodium carbonate breakthrough curves were obtained by plotting dimensionless 
exit concentration of sulphur dioxide against dimensionless time. Kinetic data obtained 
from the thermogravimetric experiments (Bhaskar, 1994) were used to simulate the 
breakthrough curves using the formulated model. 

5.2.1 Model Predicted Breakthrough Curves at Different Tempretures 

Breakthrough curves were simulated using the formulated model at different tempretures 
of 150°C, 200°C and 250°C for the sulphation of active sodium carbonate in a packed bed 
with inlet sulphur dioxide concentration and space velocity being 3500ppm and 900hr‘' 
respectively. It might be mentioned that rate constant, effective diffosivit}' and grain 
structures are the parameters which are affected by the tempreture. With increase in 
tempreture dimensionless length of the reactor, y*, and reaction modulus, cr . increases. 
The reaction modulus, cr, can be considered as a dimensionless particle size. Particle 
with small values of cr might be expected to react rapidly compared with the particle of 
large a due to less hinderence of reaction from intraparticle diffusion in smaller particle 
(Evans and Song, 1974), but the overall effect of increase in tempreture is to cause more 
removal of sulphurdioxide in the reactor as can be seen from Figure 5.1. 
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Figure 5.1: 


Pridicied Breakthrough Curves for Sulphation of Active Sodium Carbonate at 

Different Tempretures with Inlet Concentration 3 SOOppm. r.'-- / -x 
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Figured .2: Pridicted Breakthrough Curves for Sulphation of Active Sodium Carbonate at 
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5.2.2 Model Predicted Breakthrough Curves at Different Concentrations 

Breakthrough curves were simulated using the formulated model at different 
concentrations of sulphur dioxide i.e 1500, 6000 and 18000ppm for the sulphation of 
active sodium carbonate with tempreture and space velocity being 200®C and 900hf' 
respectively. With the increase in the inlet concentration of sulphur dioxide, the 
breakthrough curve is quicker as can be seen from Figure 5.2. This is also evident from 
the fact that the fractional conversion of sorbent increases with increase in the inlet 
concentration of sulphur dioxide from 1500 to 18000ppm, indicating that with increase in 
partial pressure of sulphur dioxide there is higher conversion(Bhaskar. 1994). 
Dimensionless distance along the flow, y*, and reaction modulus are not affected by 
varying the concentration at a given tempreture, however, the dimensionless time. t’. 
increases with the increase in concentration, hence breakthrough of sulphur dioxide will 
be quicker with increase in concentration. 

5.3 Validation of Simulated Breakthrough Curves 

Once the theoretical breakthrough curves for the sulphation of active sodium carbonate 
have been developed it is essential to validate them by comparing the results obtained 
from the model with those available in literature. Simulated breakthrough curves 
corroborated fairly well with the experimental results (Bhaskar, 1994) for active sodium 
carbonate and sulphur dioxide system. At tempreture 150°C and concentration 3500ppm 
it was seen from the Figure 5.3 that simulated breakthrough curve matched fairly well 
with the experimental results. For the other two tempretures (200°C and 250°C) with 
same inlet sulphur dioxide concentration it can be seen from Figure 5.4 and Figure 5.5 
that the experimental results are lying below the simulated curves indicating more 
removal in the real situation compared to the simulated one. The error between the model 
predicted and experimental breakthrough curves at concentration 3500ppm and 
tempretures 150°C, 200°C and 250°C were found out to be between 5 to 12%. 
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Dimensionless Time 


Figure o.3: Comparison or E.xperimental and Simulated Breakmrough Cur\-es for 
Suiphation of .Active Sodium Carbonate. 

* Source: Bhaskar. 1994. 





Dimensionless Time 


Fisure 5.4: Comparison of E.xperimental and Simulated Breakthrough Curves for 
Suiphation of Active Sodium Carbonate. 

* Source: Bhaskar. 1994. 
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Figure 5.5; Comparison of E.xperimentai and Simulated Breakthrough Cun.-es for 
Suiphation of .A.ctive Sodium Carbonate. 

* Source: Bhaskar. 1994. 



Figure 5.6; Comparison of E.xperimental and Simulated Breakthrough Cur\'es for 
Suiphation of Active Sodium Carbonate. 

* Source: Bhaskar, 1994. 
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Dimensionless i ime 


Figure 5.7; Companson of Experimenxai and Simuiated Breakthrough Cun.'es for 
SuiDhation of Active Sodium Carbonate. 

* Source: Bhaskar. 1994. 
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Figure 5.8: Comparison of Experimental and Simulated Breakthrough Cur^’es for 
Sulphation of Active Sodium Carbonate. 

* Source: Bhaskar, 1994. 
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It is evident that the result of graphical integration of the area above the curve to estimate 
sulphur dioxide removal based on the predicted breakthrough curv^es are on the 
conservative side for the reactor design. 

At concentration ISOOppm and tempreture 200°C it can be from Figure 5.6 that the 
simulated breakthrough curve is lying above the experimentally obtained results. 
Although the estimate for sulphur dioxide removal is on the over conservative side hut 
the reactor design becomes uneconomical for this case. At high concentrations of 
6000ppm and ISOOOppm at tempreture 200°C it can be seen from the Figure 5.7 and5.8 
that experimental results and simulated breakthrough curves are close to each other. 



6. CONCLUSIONS 


Model predicted dimensionless outlet concentration from the packed bed reactor model 
using the kinetic parameters obtained from the themogravimetric experiments performed 
by Bhaskar, 1994 corrborated fairly well with the actual dimensionless concentration for 
active sodium carbonate and sulphur dioxide system available in the litremre fBhaskar. 
1994). 

Breakthrough curves obtained by using the formulated model based on combined pore 
diffusion and chemical reaction as rate limiting steps predicted less removal of sulphur 
dioxide as compared to acmal removal at inlet sulphur dioxide concentration of 3500ppm 
and tempreture range between 200 and 250°C. It is evident that the reactor design for 
these situations based on the predicted curv'es will be on the conservative side. Infact, 
model data matched fairly well with the experimental data at 150°C. 

Model predicted breakthrough curves at different concentrations and given tempreture 
(1500, 6000 and ISOOOppm at 200°C) were found closer to the experimental results at 
higher concentrations as compared to lowest concentration. It is thus evident that the 
reactor design based on the formulated model will be on the over conserv'ative side at low 
inlet sulphurdioxide concentration. 

Errors between the predicted and experimental breakthrough curves were found to be 
between 5 and 12%, except at low concentration when design becomes over consen'ati\-e. 



7. SUGGESTIONS FOR THE FUTURE WORK 


A more comprehensive model taking into account the structural changes may be used to 
predict the conversion of sorbent with respect to time, which in turn may be incarporated 
in the packed bed reactor model to predict breakthrough curves with better accuracy. Also 
the present work did not incarporate nonisothermal conditions which may develop in the 
packed bed reactor. Work may be taken to model the packed bed reactor for predicting 
breakthrough curves under nonisothermal conditions. 
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Appendix — A 


Solution by Laplace 'IVansformation 


The original equation is 


dv di: 
dy^ ^dt- 


3 = 0 


with initial and boundary' conditions as t-fO.f*' = 1 = /(r) and fiy'.O) 
Applying Laplace Transformation. 


where. 


£ 


' dv dv 
,dy' df‘ 



= £{ 0 ] . 


^(|p) ='P-V.5)-f(s-.0) 
£(i) = 31 P 


Substituting the above values in Eq. 11). we get 


^ - aPv ^ J/P = 0 

dr 


which can be wwitten as 


dv 


dy’ 


{aP^-v -3 P ' 

The above equation can be further simplified as 

^ln(QP'’jAJ) = -^y* + C 
at y* = 0, 0(0. P) = 1/P. Substituting it in Eq. (1). we get 


1 

aP2 


ln(QP + 3) = C 


Putting this in Eq. (1) we get 


1 


\n{aP^ + ^) = —Sy' + 


1 

aP2 


ln(aP + .5) 


o = fify). 

( 1 ) 
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which can be ^litten as 

1 , f (aP^v _ 1 . . 

qP2'°\ iaP^J) j 


or. 


Therefore. 


; qP-'J - i) 

VQP - i) 


= e 


-aPy- 


(:QP-i)e-'^^-='' - 3 ] 

aP2 

g-^y- J I g-aPy- | j 

P ~ a I P2 J ~ ^ 


To derennine ihe '.-aiue of we apply inverse Lapiace Transformation on the above 
equation. Therefore 


t,- = £-'(v) 

= Uay(t’) -r -(£* - ay*)6V(i*) - -t* 

a a 


or. 



a 



-ay- 



where. Uay-{P) is a uniu' function as is given as 


Uarrin = 


0 if t’ < ay 

1 if f* > ay 



Fractional Conversion, 


Appendix B 


X 











Figure B.l: Predicted Conversion and Experimental Data for Sulphation of 
Active Sodium Carbonate. 

*Source: Bhaskar, 1994. 



Appendix C 


♦DATA USED FOR THE PREDICTION OF BREAKTHROUGH CURVES FOR SULPHATiON OF 
ACTIVE SODIUM CARBONATE. 

(♦Source: Bhaskar, 1994). 


Length of the packed bed (y) = 39.1cm. 

Radius of particle (Rq) = 0.635cm. 

Molar density of the solid reactant (/?£) = 0.0239g-mole/cm^ 


Grain shape factor (Fg )= 3. 

Particle shape factor (Fp) = 3. 

Porosity of the bed ( e,, ) = 0.42 

Superficial velocity (U) = 4cm/sec 

Fractional volume of reactive solids (w) = 1/40 




Table C.l : Physical Parameters of Active Sodium Carbonate. 





150 

0.49 

0.55 

3.14 

200 

0.47 

0.54 

2.70 


250 


0.45 


0.53 


2.50 
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Table C.2: Values of Rate Constant (k) and Effective Diffusivit}' (Dj) for Sulphation of 
Active Sodium Carbonate. 







Appendix D 


Table D.l: Experimental and Pridicted Values of Breakthrough Curves for Sulphation of 
Active Sodium Carbonate at 150®C with Inlet Sulphurdioxide Concentration of 
3500ppm. 


Dimensionless Outlet Concentration 

Dimensionless Time Experimental* Predicted 


0.000 

0.000 

0.000 

0.011 

0.462 

0.627 

0.022 

0.678 

0.709 

0.033 

0.791 

0.754 

0.044 

0.820 

0.783 

0.057 

0.873 

0.804 

0.067 

0.888 

0.821 

0.083 

0.876 

0.844 

0.100 

0.812 

0.853 

0.117 

0.885 

0.865 

0.150 

0.885 

0.883 

0.167 

0.888 

0.887 

0.184 

0.906 

0.895 

0.200 

0.864 

0.901 


♦Source: Bhaskar, 1994. 
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Table D.2: Experimental and Pridicted Values Breakthrough Curv^es for Suiphation of 

Active Sodium Carbonate at 200°C and Inlet Sulphurdioxide Concentration of 
3500ppm. 


Dimensionless Time 


Dimensionless Outlet Concentration 
Experimental" predicted 


0.000 

0.000 

0.000 1 

0.017 

0.426 

0.528 I 

0.034 

0.626 

0.665 

0.051 

0.683 

0.721 

0.068 

0.747 

0.757 

0.085 

0.753 

0.782 

0.102 

0.730 

0.801 

0.127 

0.753 

0.823 

0.153 

0.761 

0.838 

0.178 

0.822 

0.852 

0.208 

0.788 

0.861 

1 

0.229 

0.788 

0.870 

0.255 

0.805 

0.877 

0.280 

0.813 

0.886 


^Source: Bhaskar, 1994. 
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Table D.3: Experimental and Predicted Values Breakthorugh Curves for Sulphation of 

Active Sodium Carbonate at 250°C and Inlet Sulphurdioxide Concentration of 
SSOOppm. 


Dimensionless Outlet Concentration 
Dimensionless Time Experimental * Predicted 


0.000 

0.000 

0.000 

0.025 

0.153 

0.376 

0.051 

0.361 

0.541 

0.076 

0.465 

0.620 

0.101 

0.572 

0.670 

0.126 

0.634 

0.706 

0.152 

0.685 

0.732 

0.189 

0.712 

0.761 

0.227 

0.685 

0.784 

0.265 

0.715 

0.800 

0.303 

0.754 

0.815 

0.341 

0.768 

0.827 

0.379 

0.805 

0.837 

0.417 

0.826 

0.848 


♦Source; Bhaskar, 1994. 




Table D.4: Experimental and Predicted Values of Breakthrough Curv'es for Sulphation of 
Active Sodium Carbonate with Inlet Sulphurdioxide Concentration of 
ISOOppm and tempreture 200°C. 

Dimensionless Outlet Concentration 

Dimensionless time Experimental* Predicted 


0.000 

0.000 

0.000 

0.014 

0.233 

0.525 

0.021 

0.358 

0.591 

0.028 

0.417 

0.636 

0.035 

0.417 

0.669 

0.042 

0.512 

0.695 

0.055 

0.501 

0.733 

0.065 

0.565 

0.747 

0.076 

0.611 

0.771 

0.087 

0.678 

0.780 

0.098 

0.717 

0.797 

0.110 

0.717 

0.810 

0.120 

0.724 

0.816 

0.131 

0.748 

0.824 


* Source bhaskar, 1994. 
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Table D.5: Experimental and Predicted Values of Breakthrough Curv’es for Sulphation of 
Active Sodium Carbonate with Inlet Sulphurdioxide Concentration of 
bOOOppm and Tempreture 200°C. 

Dimensionless Outlet Concentration 

Dimensionless Time Experimental* Predicied 
















































* Source; Bhaskar, 1994. 
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Table D.6: Experimental and Predicted Values of Breakthrough Curves for Sulphation of 
Active Sodium Carbonate with Inlet Sulphurdioxide Concentration of 
ISOOOppm and tempreture 250°C. 

Dimensionless Outlet Concentration 

Dimensionless Time Experimental* Predicted 


0.000 

0.000 

0.000 

0.087 

0.863 

0.803 

0.174 

0.910 

0.857 

0.261 

0.958 

0.882 

0.348 

0.952 

0.899 

0.435 

0.943 

0.911 

0.522 

0.964 

0.920 

0.654 

0.960 

0.929 

0.783 

0.976 

0.937 

0.915 

0.940 

0.945 

1.044 

0.990 

0.953 

1.178 

0.965 

0.961 

1.305 

0.971 

0.971 

1.570 

0.982 

0.979 


/ 


"^'Source: BhaskarJ994. 




Appendix E 

SOLUTION FOR THE FRACTIONAL CONVERSION, USING RUNGE- KUTTE METHOD 

implicit real (a-z.o-z) 

print*, ’specify sigma,h,tmax’ 

read* ,sigma,h,tmax 

open(unit=2 1 ,riIe=concen.out’) 

ti=0.0 

x=0.0 

n=int(tmax/h) 
write(2 1 ,24)sigma 

24 format( 5x,'the value of sigma =‘,2x,f6.3/) 

WTile(21,24) 

25 format(8x,'time t*,7x,’beta’,9x,’x(t*)’y) 
xi=0.0 

bit=beta(xi,sigma) 

t=0.0 

write(21,22)y,bit,xi 

do i=2,n 

xi=x 

rk.O=h*beta(xi,sigma) 

rk 1 =h* beta(xi+0.5 * rkO.sigma) 

rk2=h*beia(xi+0.5*rkl, sigma) 

rk3=h*beta(xi+rk2,sigma) 

x=xi+(rk0+2.0*rkl+2.0*rk2+rk3)/6.0 

b=beta{x, sigma) 

bit=b 

t=h*(i-l) 

if(x.ge.0.999) goto 90 
write(2 1 ,22)t,bit,x 
22 format(5x,4(fl0.6,3x)) 

100 end do 
90 write(21,22)t,bit,x 

stop 
end 

function beta(x, sigma) 
implicit real(a-h,o-z) 
exp— 1/3 
expo2— 2/3 

fl =2.0* sigma* sigma* (( 1 .0-x)* *expo 1 - 1 .0) 

f2=((1.0-x)**expo2)/3.0 

beta=1.0/(n+f2) 

return 

end 



